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ABSTRACT 

We discuss free- free radio emission from ionized gas in the intergalactic medium. Because the emissivity 
is proportional to the square of the electron density, the mean background is strongly sensitive to 
the spatial clumping of free electrons. Using several existing models for the clumping of ionized gas, 
we find that the expected free-free distortion to the cosmic microwave background (CMB) blackbody 
spectrum is at a level detectable with upcoming experiments such as the Absolute Radiometer for 
Cosmology, Astrophysics, and Diffuse Emission (ARCADE). However, the dominant contribution to the 
distortion comes from clumpy gas at z ^ 3, and the integrated signal does not strongly constrain the 
epoch of reionization. In addition to the mean emission, we consider spatial fluctuations in the free-free 
background and the extent to which these anisotropics confuse the search for fluctuations in 21 cm line 
emission from neutral hydrogen during and prior to reionization. This background is smooth in frequency 
space and hence can be removed through frequency differencing, but only so long as the 21 cm signal 
and the free-free emission are uncorrelated. We show that, because the free-free background is generated 
primarily at low redshifts, the cross-correlation between the two fields is smaller than a few percent. 
Thus, multifrequcncy cleaning should be an effective way to eliminate the free-free confusion. 

Subject headings: cosmic microwave background — large scale structure of Universe — diffuse radiation 



1. INTRODUCTION 

Plans for upcoming low-frequency radio experiments 
able to measure the 21 cm background associated with 
neutral hydrogen during, and prior to, reionization era 
(e.g., Scott & Rees 1990; Madau et al. 1997; Zaldar- 
riaga et al 2003; Morales & Hewitt 2003) have motivated 
study of the foregrounds that may contaminate such mea- 
surements. Synchrotron emission from the Milky Way, 
low frequency radio sources (Di Matteo et al 2002), and 
free-free emission from free electrons in the intergalactic 
medium (Oh & Mack 2003) are now thought to be the 
chief sources of confusion. 

While the free- free background contaminates 21 cm 
studies, the emission itself captures important physics of 
the ionized component of the intergalactic medium (IGM) . 
In particular, because the emissivity is proportional to 
the square of the electron number density, free-free emis- 
sion is strongly sensitive to whether electrons are spatially 
clumped or distributed smoothly in the IGM. Initial esti- 
mates of the free-free background suggest that gas clump- 
ing boosts the specific intensity by a factor from order 
unity at ~ 20 to over 100 at z ^3 (Oh 1999). 

Here, we consider the direct detection of free- free emis- 
sion through the distortion it creates in the cosmic mi- 
crowave background (CMB) blackbody spectrum at low 
radio frequencies (Loeb 1996; Oh 1999). We calculate the 
mean distortion temperature using a variety of both ana- 
lytic and numeric clumping models already existing in the 
literature. While the expected brightness varies by factors 
of a few between the different models, we suggest that 
the distortion to the CMB is at a level detectable with 
the planned Absolute Radiometer for Cosmology, Astro- 
physics, and Diffuse Emission (ARCADE; Kogut 20003) 
experiment. A detection of the distortion could help to 
discriminate between existing clumping models and to con- 
strain the integrated ionization history of the universe. We 



also show that ionized halos at z 3 dominate the back- 
ground, with only a relatively insignificant contribution 
from the reionization era. 

In addition to the mean background, we estimate the 
magnitude of spatial fluctuations in the free-free intensity 
that appear because the ionized clumps are a biased tracer 
of the dark matter density field. We find that the fluctu- 
ations are comparable to or greater than those from the 
21 cm signal on the relevant scales. Fortunately, as shown 
by Zaldarriaga et al (2003), the smoothness of the free- 
free background in frequency space should allow one to 
clean these sources in 21 cm maps (see also Morales & 
Hewitt 2003); the cleaning is quite efficient but relies on 
the foreground and 21 cm signals being uncorrelated. In 
fact the two should be anticorrelated, because the free-free 
emission comes from ionized halos while the 21 cm signal 
comes from neutral regions. Here we show that, because 
nearly all of the free- free background arises at z ^ 3, the 
anti-correlation is at the level of few percent, suggesting 
that adequate cleaning is possible. 

The discussion is organized as follows: in §2, we brieffy 
discuss free-free emission from ionized halos and its po- 
tential detectability as a distortion to the CMB. In §3, we 
discuss the detectability of spatial fluctuations in the free- 
free background and the extent to which free-free emission 
contaminates 21 cm measurements. Throughout the pa- 
per, we make use of the WMAP-favored LCDM cosmolog- 
ical model (Spergel et al 2003). 

2. THE MEAN FREE-FREE SIGNAL 
First, we write the free-free emission coefficient as 



= 5.4 X lO-^^nlT-'/^gff{iy,T,)e-''''/ 



ergs cm ^ s ^ 



Hz"^ sr^^ , 



(1) 



where the Gaunt-factor can be approximated in the radio 



regime as gff{v,Tf.) ~ 11.96Tj 
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(Lang 1999). The 
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Fig. 1. — (a) The clumping factor of electrons in the IGM as 
a function of redshift. These clumping factors come from Haiman 
et al (2001) (long-dashed line: a reionization model that includes 
minihalos as well as star- forming halos with Tv > 10* K), Ben- 
son et al (2001) (dashed line: semi-analytic model including cool 
neutral gas, or galaxies, in ionized halos as well as the clumping 
of gas outside halos), and Gnedin & Ostriker (1997) (dot-dashed 
line: a direct measurement in numerical simulations). The dot- 
ted line is the clumping factor presented in Oh (1999), while the 
solid line is a similar calculation described in the text based on 
the star formation history of the universe, (b) The mean bright- 
ness temperature of free-free emission at an observed frequency of 
2 GHz. Here, we show both the differential (thin lines) and cu- 
mulative brightness temperatures (thick lines) as a function of the 
redshift. The curves refer to the clumping models in the top panel. 
For reference, we show the case where C{z) = 1 and the ionized 
gas distribution is taken to be spatially smooth. We also show the 
current limit on the free-free distortion of the CMB (top line with 
Yff = T^/TcMsihu/kTcMB)'^ < 1-9 x 10"^ K; BersaneUi et al 1996) 
and the expected constraint from the ARCADE experiment (Kogut 
2003). 

cumulative specific intensity is simply = / dxiv{z)x1{z) 
C{z)/{1 + z)'^, where x is the conformal time or lookback 
distance from the observer and = fobs(l + z) with 
Vohs the observed frequency. The electron fraction Xe{z) 
captures the ionization history of the IGM. In our calcula- 
tions, unless specified as part of the clumping description, 
we take a model for reionization based on the UV radiation 
from star formation (e.g., Haiman & Holder 2003; Chen 
et al 2003), with a cumulative optical depth to electron 
scattering of ~ 0.14 (consistent with the recent measure- 
ments from CMB polarization data; Kogut et al 2003). We 
will quote results in terms of the brightness temperature 
Tfc — c^Iy /2v1^^kB (vahd for the Rayleigh- Jeans part of 
the spectrum). 

In Eq. (1), we have introduced C{z) — {nl)/{ne)^, the 
clumping factor of the electron density field. As has been 
noted previously (Oh 1999), the free-free intensity depends 
strongly on this factor. Instead of choosing a specific pre- 
scription, we make use of a variety of published models for 
C{z). One approach is to estimate the clumping from the 
large-scale matter distribution (by, for example, assum- 



ing that ionized electrons reside in virialized dark matter 
halos). Several such models are shown in Fig. 1(a), in- 
cluding estimates based on analytic (Haiman et al 2001), 
semi-analytic (Benson et al 2001), and numerical (Gnedin 
& Ostriker 1997) techniques. These clumping models differ 
in detail, especially around z ~ 10, because they include 
different aspects of the reionization process, such as the 
presence of partly ionized minihalos (Haiman et al 2001). 

The second approach, taken by Oh (1999), is to con- 
struct a model for the production rate of ionizing pho- 
tons TT-ion and assume that, throughout the universe, ion- 
ization equilibrium is a good approximation. The quan- 
tity XgC(z) is then fixed through the relation njon = 
Xg(z)rig(z)C(2:)aroc, where arcc is the recombination co- 
efficient. We show the clumping factor estimated by Oh 
(1999) in Fig. 1(a) as a dotted line; at high redshifts it 
assumes that the ionizing photon production rate is pro- 
portional to the mass in halos with virial temperatures 
> 10^ K (with ~ 400 ionizing photons per collapsed 
baryon, appropriate for a star formation efficiency of 10% 
and Population II stars with a standard IMF). We also 
show a similar estimate using the analytical model of Hern- 
quist & Springel (2003) for the star formation history of 
the universe, again assuming Population II stars with a 
standard IMF. While the star formation history is model 
dependent at high redshifts, this model fits observational 
measurements at z 3 (e.g., Somervillc, Primack, & Faber 
2001) reasonably well. The Oh (1999) approach predicts 
stronger clumping than models based on models of the 
IGM gas do, suggesting that the latter approach may un- 
derestimate the dumpiness of galaxies. However, all of the 
clumping models generally have similar evolutionary his- 
tories, with C{z — 20) ~ 1 increasing to C(2; = 0) ~ lO'^. 
We will find below that the detailed behavior at z ^ 10 
has little effect on the observable results. 

As shown in Fig. 1(b), dT^/dz closely traces the red- 
shift evolution of C{z). In a smooth IGM, with C{z) = 1, 
a wide range of redshifts contribute approximately equally 
to the free-free background, decreasing only around reion- 
ization (z ~ 15 in our model). Clumping increases the 
contribution to the free-free background by a factor C{z) 
Sit each redshift: because of the sharp increase in the 
clumping factor at z 3, most of the free-free background 
arises at these low redshifts and not during the reioniza- 
tion epoch. Unfortunately, since one measures only the in- 
tegrated background, the redshift evolution of the clump- 
ing factor is not directly observable. Note that the mod- 
els considered give predictions that vary by a factor of a 
few. This suggests that these measurements may provide 
a discriminant between some of these models, especially 
between the two different approaches to clumping calcula- 
tions. We also see that, independent of the detailed mod- 
eling, the predicted distortion to the CMB spectrum is 
at a level that can be detected with ongoing experiments, 
such as ARCADE (Kogut 2003). While such a detection 
would not allow us to extract with confidence a specific 
model for C (z) or to constrain the reionization history in- 
dependently, a detection or upper limit of ^ 10^"^ K at a 
frequency of 2 GHz would place useful limits on current 
models of the ionized gas distribution. 

In these calculations, we have assumed Tg = 2 x 10"* 
K, consistent with measurements of the Lya forest at z 
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2.4-3.9 (Zaldarriaga et al 2001). The free-free emission 
does not depend strongly on T^, {e„ oc T^°-^^, including 
the Gaunt factor), so allowing the temperature to vary 
as a function of density, halo mass, or redshift only leads 
to a reduction in the mean brightness by a few tens of 
percent or less. Thus, we do not expect uncertainties in 
the thermal history of the ionized gas to be a significant 
limitation in these estimates. 

3. ANGULAR FLUCTUATIONS FROM FREE-FREE 
EMISSION 

We now consider the observability of spatial fluctuations 
in the free-free intensity. In general, fluctuations in the 
brightness temperature of the free-free background arise 
from fluctuations in the electron density field. For sim- 
plicity, we will neglect anisotropics generated by fluctu- 
ations in the ionized fraction during reionization, as the 
fractional contribution to the mean intensity is dominated 
by ionized gas at z 3. 

We define the angular power spectrum of the free-free 
background to be {Tb{l)Tb{l')) = {27Tf5{l + l')Ci, , where 
Tii{l) is the Fourier transform of free- free brightness tem- 
perature. Using the Limber approximation (Limber 1954), 
we can write the brightness temperature anisotropy power 
spectrum, observed at a frequency v, as 



cr = 



dx 

"■A 



(1 + Z)4 



n 2 



P.. 
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where (Ia is the comoving angular diameter distance and 
Pee (k) is the three-dimensional power spectrum of the elec- 
tron field. To calculate this quantity we make use of the 
halo approach to large scale structure clustering (Cooray 
& Sheth 2002). We assume that on the scales of interest 
(from degrees to arcminutes), we can describe fluctuations 
in the density field through the 2-halo part of the power 
spectrum, which includes the clustering of different ion- 
ized halos. Thus we write Pee(k) = b^^P^^'^{k) in terms of 
the linear power spectrum. We assume that the ionized 
patches exist in all halos with a minimum mass, Mj^iniz) 
that corresponds to T„ = 10** K and calculate the bias fac- 
tor of these halos with respect to the linear density field. 
Thus 

APbi,,,,{M,z)dn/dM 
lM,^,„M^dn/dM ' 

where 6haio(-^^, z) is the bias factor of the dark matter ha- 
los (Mo et al. 1997). Note the additional IVP weighting, 
which takes into account the fact that the relevant bias 
factor is that of ti^. Essentially, we assume that the free 
electrons are in the form of ionized clumps in each dark 
matter halo, and we assign a halo occupation number pro- 
portional to the dark matter halo mass (in analogy to the 
halo approach for galaxy biasing). With this weighting, 
bee exceeds the dark matter bias, although the difference 
is minor at redshifts of a few. 

Our results for the angular fluctuations are summarized 
in Fig. 2(a). The typical fluctuations are 10^^ K at 2 
GHz, a few percent of the mean brightness. This is consis- 
tent with a mildly biased field that traces the large scale 
dark matter field and has its dominant contribution at red- 
shifts of a few. Again the Oh (1999) model predicts sub- 
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Fig. 2. — (a) The rms brightness temperature fluctuations 
of the free-free background at !/ot,s = 2 GHz, where ATrms = 
^1{1 + 1)/2tvCi. The curves refer to the clumping models in Fig. 1. 
The thin lines are for fluctuation power based on equation (2) using 
the bias factor in equation (3) and the linear power spectrum, while 
the thick lines include non-linearities in the dark matter density fleld 
following Smith et al (2003). At Z > 10^, non-linearities are impor- 
tant and our estimate for fluctuation power spectrum is likely to 
be only approximate given uncertainties in the spatial distribution 
of ionized clumps within each halo. At these small scales, we may 
also have underestimated power as we do not consider the Poisson 
contribution related to the flnite number density of free- free sources. 
For reference, we also show primordial CMB anisotropics as well as 
noise power spectra for ARCADE and for the SKA, assuming one 
month of continuous observations for the latter, (b) The brightness 
temperature fluctuations at u^hs = 140 MHz. Here, we show the 21 
cm signal in a narrow window around z = 10 (taken from Zaldar- 
riaga et al 2003), the CMB, and free-free fluctuations. The dotted 
line shows the free-free fluctuations including only those halos in the 
redshift interval of the 21 cm power spectrum; the long-dashed line 
shows the (absolute value of the) cross-correlation between the 21 
cm and free-free signals in this redshift window. 



stantially larger fluctuations than the other models, sug- 
gesting that measurements may distinguish between the 
two predictions (though without constraining the details 
of the large scale IGM clumping) . For comparison, we also 
show the noise level, at each multipole, expected for the 
ARCADE experiment, if it attempts to make anisotropy 
measurements with a noise per beam size pixel of 100 fiK 
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and a beam size of 16°, and for the planned Square Kilome- 
ter Array (SKA). In calculating the SKA noise, we follow 
the approach of Zaldarriaga et al (2003) and assume a sys- 
tem temperature Tgys = 50 K at 2 GHz, a bandwidth of 10 
MHz, and continuous observations over 4 weeks. We also 
show the anisotropy power spectrum of primordial CMB 
fluctuations. While the CMB fluctuations swamp those 
expected from the free-free background, the CMB contam- 
ination can be reduced through two methods: (1) because 
free-free fluctuations are dominated by point-like clumped 
regions in the IGM, which are likely to be associated with 
dark matter halos, one can fllter the large scale fluctua- 
tions of CMB and concentrate on the remaining point-like 
emission, or (2) one can use existing cosmic-variance lim- 
ited CMB maps at higher frequencies to clean the CMB 
contamination. We see that SKA can in principle detect 
free-free fluctuations on scales I ^ 3000. The main source 
of confusion will be the uncertain contribution from extra- 
galactic radio sources, such as radio galaxies (Di Matteo et 
al 2002). Because these soiirces arc also point-like, separat- 
ing them from free-free fluctuations will require frequency 
information. 

To estimate the background for 21 cm measurements, we 
consider free-free fluctuations at z^obs = 140 MHz, corre- 
sponding to Hne emission from neutral hydrogen at z ss 10. 
We summarize our results in Fig. 2(b). We also show the 
SKA noise power spectrum, where we have set Tgys = 200 
K and the bandwidth to 0.2 MHz, because one wants nar- 
row redshift slices in order to observe many independent 
regions for 21 cm fluctuations. The free-free fluctuations 
are comparable to or greater than the rms fluctuations in 
the 21 cm background, as estimated by Zaldarriaga et al 
(2003), while contamination by the CMB is negligible. As 
noted by Oh & Mack (2003), free- free fluctuations present 
a substantial foreground; however, the free-free spectrum 
is smooth in frequency space, while the 21 cm signal varies 
strongly over ~ 0.2 MHz. Thus multi- frequency differenc- 
ing can efficiently remove the smooth background (Zaldar- 
riaga et al 2003), but only so long as the two signals are 
uncorrelated. We expect this to be the case with free-free 
emission: halos with ionized gas will emit free-free radia- 
tion but not 21 cm radiation, so the two fields are strongly 
anti-correlated. 

To test how important this effect is, we show the free- 
free emission from halos in the same redshift slice as the 
21 cm signal (a window of 0.2 MHz about z = 10) with 
the dotted line. To construct this curve we have used the 
Oh (1999) clumping model; this gives a mean ionizing rate 
roughly similar to the simulations on which the simple es- 
timate of Zaldarriaga et al (2003) are based (Sokasian et 
al 2003). Note that AT^ms oc xlC(z) and can thus be 
easily scaled for different clumping scenarios. In partic- 
ular, in the Oh (1999) approach the free-free signal will 
just be proportional to nion- We then estimate the cross- 
power spectrum between the free- free and 21 cm signals by 
assuming that they are perfectly anti-correlated on large 
scales; the (absolute value of the) result is shown by the 
long-dashed line in the Figure. This curve is thus our esti- 
mate for the maximally cleaned 21 cm map; it shows that 
the two are correlated only on the percent level. Because 
it is smaller than the noise power spectrum for SKA ob- 
servations, the free- free signal can be efficiently removed 
with multi-frequency data. 



4. CONCLUSIONS 

To summarize, we find that the cumulative specific in- 
tensity of the free-free background at low radio frequen- 
cies is strongly sensitive to the spatial clumping of free 
electrons in the IGM. Using a variety of existing models 
for the clumping of ionized gas and its redshift evolution, 
we find that the expected free-free distortion to the CMB 
blackbody spectrum, at frequencies of order a few GHz 
and below, is within the reach of upcoming experiments 
such as ARCADE. The dominant contribution to the back- 
ground is from ionized gas at 2; ~ 3, because the clumping 
factor is generally large below this redshift. The free- free 
background varies across the sky at the level of a few to 
at most ten percent. These fluctuations can be detected 
with experiments such as the SKA, though a careful sepa- 
ration of CMB anisotropics and other low-frequency radio 
point sources will be required. At frequencies of order 200 
MHz, where observations of 21 cm radiation from neu- 
tral hydrogen during and before the epoch of reionization 
are planned, angular fluctuations from free-free emission 
could become a significant source of confusion. However, 
we have shown that the cross-correlation between the 21 
cm and free-free signals is small, suggesting that the free- 
free foreground can be removed to good accuracy. 
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